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We study few-body problems in mixed dimensions with N > 2 heavy atoms trapped individuahy 
in parahel one-dimensional tubes or two-dimensional disks, and a single light atom travels freely in 
three dimensions. By using the Born-Oppenheimer approximation, we find three- and four-body 
bound states for a broad region of heavy-light atom scattering length combinations. Specifically, 
the existence of trimer and tetramer states persist to negative scattering lengths regime, where no 
two-body bound state is present. These few-body bound states are analogous to the Efimov states 
in three dimensions, but are stable against three-body recombination due to geometric separation. 
In addition, we find that the binding energy of the ground trimer and tetramer state reaches its 
maximum value when the scattering lengths are comparable to the separation between the low- 
dimensional traps. This resonant behavior is a unique feature for the few-body bound states in 
mixed dimensions. 

PACS numbers: 34.50-s, 03.65.Ge, 34.50.Cx 



I. INTRODUCTION 

One striking feature of few-body physics is the presence 
of universality under a resonant short-range interaction, 
where the low-energy behavior of the system does not 
depend on the details of its structure or interactions at 
short distances. Of particular interest is the existence of 
bound trimer states for three identical bosons in three 
dimensions with a resonant two-body interaction, as dis- 
cussed in 1970 by Vitaly Efimov [1]. At infinite scattering 
length, these three-body bound states form an infinite 
geometric spectrum with a constant ratio between two 
successive binding energies, indicating a discrete scaling 
symmetry Besides, the bound trimer states persist 
rather counterintuitively to negative scattering length 
regime, where two-body bound states are not existent. 
After its original proposal, the Efimov physics has at- 
tracted great attention in multi-disciplinary systems, in- 
cluding atomic nuclei d, ^He trimer s [H, Q, and other 
molecules 0- However, a direct evidence of such peculiar 
behavior was not achieved for more than three decades 
until its first observation in an ultracold gas of neutral 
atoms [H. Thanks to the extraordinary controllability 
of the mutual atomic interaction by tuning through a 
magnetic Feshbach resonance, signatures of trimer bound 
states have been observed in trapped atomic gases for 
both negative and positive scattering length regimes [sl- 

la. 

In addition to the original problem of identical bosons, 
the study of three-body physics has been extended to 
a variety of other three-particle systems [l3|, includ- 
ing three distinguishable particles with different scatter- 
ing lengths and/or different masses ^14nl6j, two iden- 
tical fermions with a third atom [l5|, [l8|, and the 
three-atom systems with with non-zero angular momen- 
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tum [l9|. Of particular interest is the case of three 
distinguishable fermions in an ultracold gas of three- 
component ^Li atoms. In such a system, there exits a 
broad magnetic Feshbach resonance such that all three 
scattering lengths can be tuned around resonance simul- 
taneously [20j, leading to a promising candidate to ob- 
serve the few-body universal behavior [21-29]. Besides, 
the few-body problem has also been analyzed in ultracold 
gases of different atomic species by tuning the interaction 
across an interspecies Feshbach resonance (3Q| . 

Due to the multi-channel nature of the inter-atomic 
interaction, the Efimov states in the three-dimensional 
(3D) ultracold atomic gases are only metal-stable states. 
Through the three-body recombination process, two of 
the three atoms in an Efimov trimer can further form a 
deeply bound dimer and the third one would escape from 
the trap. In order to prepare stable trimer states, one has 
to figure out a mechanism to significantly reduce or even 
prevent three-body recombination. Since the three-body 
recombination process only occurs when three atoms all 
come to a close range, one possible route towards this goal 
is to use geometric confinement to separate atoms such 
that they can not travel to a same spot. For instance, if 
two of the three atoms are individually trapped in two 
spatially separated one-dimensional (ID) tubes or two- 
dimensional (2D) disks, and interact with each other via 
the third atom which is free in all three dimensions (3D), 
the three-body recombination is inherently forbidden and 
the trimer states are stable if they exist in this mixed 
dimensional configuration. 

The few-body problem in mixed dimensions has been 
recently discussed by Nishida and Tan [31], where they 
consider two species of atoms confined in different dimen- 
sions and find trimer bound states for a certain range of 
mass ratio. However, since the atoms in lower dimensions 
are not geometrically separated, this configuration suf- 
fers the same problem of three-body recombination and 
the trimer states are unstable. Therefore, Nishida con- 
sider the problem of two atoms trapped in two separated 
ID tubes or 2D layers, and interacts with the third atom 
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which is free in 3D [32|,|36j. This 1D-1D-3D or 2D-2D-3D 
mixture thus can support stable Efimov trimer states. 

In this manuscript, we adopt another approach based 
on Born-Oppenheimer approximation (BOA) to study 
some few-body problems in mixed dimensions, and in- 
vestigate the existence and properties of stable few-body 
bound states in a variety of configurations. For the three- 
body problems, we consider the systems with two heavy 
atoms trapped in two parallel ID tubes (1D-1D-3D) or 
2D disks (2D-2D-3D), plus one light atom moving freely 
in 3D (see Fig. 1 for illustration). We conclude that the 
light atom can induce an effective interaction between 
the two heavy atoms which are spatially separated by 
the low dimensional traps. Due to this effective interac- 
tion, the two heavy atoms can be bound with each other 
and lead to the formation of a three-body bound state 
in a very broad parameter region, including the regimes 
with negative 5- wave scattering lengths between the light 
and the two heavy atoms, where two-body bound states 
are not present. 

In addition to their existence in mixed dimensions, 
the universal three-body bound states also acquire some 
unique features due to the geometric confinement. Espe- 
cially, we find that the two heavy atoms experience the 
strongest effective interaction when the scattering length 
between heavy and light atoms equals to the distance be- 
tween the two low-dimensional traps. As a consequence 
of this resonance phenomenon, the binding energy of 
the ground trimer state takes a peak value around the 
resonance point, where the scattering length is of a fi- 
nite value. We emphasis that, the BOA provides a very 
clear physical picture with which the new resonance phe- 
nomenon in the mixed-dimensional systems can be easily 
explored and clearly described. 

We also compare our results with the exact expres- 
sion [32I, [sil given by an effective field theory, and con- 
clude that BOA works well even in systems with a mass 
ratio only about 6. This finding suggests that BOA is a 
powerful tool for the study of stable heteronuclear few- 
body bound states in mixed dimensions. 

To demonstrate the usage of BOA for general few-body 
problems, we consider as an example the 1D-1D-1D-3D 
geometry with three heavy atoms confined individually 
in parallel ID tubes and a light atom in 3D free space. 
We find four-body bound states living in a wide range of 
scattering lengths. A similar resonance phenomenon is 
also observed when the scattering length becomes close 
to the mutual distances between ID tubes, in which case 
the binding energy of the ground tetramer state reaches 
its maximum when the three ID tubes form an equilateral 
triangle. We also show the scheme to generalize the BOA 
to problems with > 3 heavy atoms and a single light 
one in an arbitrary mixed dimensional geometry. The 
mixed dimensional systems discussed in this manuscript 
can be realized in a mixture of two-species ultracold gases 
with species-selective dipole potential, as illustrated in 
recent experiments [33| . 

The remainder of this manuscript is organized as fol- 



lows. In Sec. II we first consider the 1D-1D-3D geom- 
etry and outline the BOA approach for the three-body 
problem. We calculate the effective interaction potential 
between the two heavy atoms, and observe the new res- 
onance phenomenon. In Sec. Ill we solve for the three- 
body bound states, from which we conclude that a sta- 
ble trimer state can exist in a broad parameter region, 
and the binding energy of the ground trimer state takes 
largest value under the new resonance condition. Similar 
results of the 2D-2D-3D system are shown in Sec. IV. 
In Sec. V we extend the usage of BOA to the four-body 
problem in 1D-1D-1D-3D geometry, and discuss the ex- 
istence and properties of bound tetramer states. In Sec. 
VI, we show the general scheme to apply BOA in prob- 
lems with more than 3 atoms in arbitrary mixed dimen- 
sional geometries. Our main findings are concluded in 
Sec. VII, and the Bethe-Peierls boundary condition used 
in our BOA approach is derived in Appendix A. 



II. BOA FOR THREE- ATOM BOUND STATES 
IN 1D-1D-3D SYSTEMS 

In this section we present the Born-Oppenheimer ap- 
proach for a three-body system with two heavy atoms 
individually trapped in two parallel ID tubes and a light 
atom moving freely in the 3D space. The straightfor- 
ward generalization to 2D-2D-3D systems will be given 
in Sec. IV, while the discussion for four-body problems 
in 1D-1D-1D-3D systems is given in Sec. V. 



A. System and Hamiltonian 

As shown in Fig. 1(a), the 1D-1D-3D system includes 
two heavy atoms Ai and plus a light atom B. The 
atoms Ai and A2 are trapped in two parallel ID tubes 
centered along the lines (x = ±L/2, = 0), while the 
light atom B moves freely in the 3D space. The quan- 
tum state of this system can be described by the wave 
function ^(r^; ^1, ^2), where zi^2 are the z-coordinate of 
atoms Ai^2 in the ID tubes, and tb = {xb^Vb^ zb) is the 
coordinate of atom B in 3D. 

In this manuscript, we use the natural units h = ttib = 
L = 1, where is the mass of atom B. The Hamilto- 
nian for the motion of the three atoms is 

2 2mi dzf 2m2 

where mi ^2 are the masses of atoms Ai^2 in the natural 
unit, and Vib,2B are the interaction potentials between 
Ai^2 and 5, respectively. In this work we only consider 
the cases where the distance L between the two tubes is 
much larger than the characteristic length of the inter- 
action potential between Ai and A2. Hence the A1-A2 
interaction can be safely ignored. 
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FIG. 1: (color online) (a) The 1D-1D-3D system with two 
heavy atoms Ai and A2 confined in two ID tubes and the 
light atom B moving freely in 3D. (b) The 2D-2D-3D system 
with two heavy atoms Ai and A2 confined in two 2D planes 
and the light atom B moving freely in the 3D space. 



B. BOA for three-body bound states 

The three-body bound state is given by the solution of 
the eigen-equation 



(2) 



When the masses of the heavy atoms ^41^2 is much larger 
than the one of B^ or mi, 2 ^ 1 in the natural unit, the 
eigen-equation (j2j) can be solved with BOA. This approx- 
imation is applicable when the motion of the heavy atoms 
Ai^2 is slow enough such that the quantum transitions 
between different instantaneous eigen-states of the light 
atom B with fixed positions Zi^2 of ^1,2 are negligible. 
Therefore, the total wave function ^ of the three-body 
bound state can be approximated as a factorized form 



(3) 



where il^ivB-, ^1, ^2) is an instantaneous bound-state solu- 
tion of the eigen-equation of the Hamiltonian of atom B 
with fixed values of zi and ^2. 

As shown in Appendix A, we can further replace the 
interaction potentials Vib and V2B with the Bethe-Peierls 
boundary conditions 

^{riB ^ 0) a f 1 - ^\ + (9(riB); (4) 

^{r2B ^ 0) (X f 1 - ^\ + 0{r2B)- (5) 

Here riB,2B are the relative distances between the heavy 
atoms Ai^2 and the light atom 5, ai,2 are the mixed- 
dimensional scattering lengths between Ai^2 and B. No- 
tice that the Bethe-Peierls boundary conditions (HI and 
(j5j) are derived from the first-principle calculation where 
the 3D motion of all the three atoms Ai^2 and B are 
taken into account. Then the mixed-dimensional scat- 
tering lengths ai,2 are determined by both the 3D s- 
wave scattering lengths between Ai^2 and 5, as well as 
the intensity of the transverse confinements of the ID 



traps. Thus, ai,2 can be tuned either through a 3D mag- 
netic Feshbach resonance [35| or via a mixed-dimensional 
confinement- induced resonance (sij . 

With the Bethe-Peierls boundary conditions, the wave 
function ip{rB^ Z2) is determined by 



-V|?A(fB, 2^1,^2) 



Kff(^l, ^2)V^(rB, 2^1, Z2), (6) 



through which the shape of the wave function 
V^lre, ^1, ^2) and the relevant eigen-energy leff(^i, ^2) can 
be determined for a given value of ^1,2- 

In the approach of BOA, the instantaneous energy 
Vef^{zi^Z2) of the light atom B serves as an effective 
potential between the two slowly moving heavy atoms. 
Then the wave function ^(^1,^2) in Eq. ([3]) satisfies the 
Schrodinger equation 



1 



1 



2mi dzf 2m2 dz2 



Vei^{zi,Z2) 



(^1,^2) 



^(/)(^1,^2), (7) 



where E is the total energy of the trimer state defined 
in Eq. (j2j). In this manuscript, we focus only on the 
ground state of the three-body eigen-equation ([2]), which 
is consisted of the ground-state solutions i/j and (j) of ^ 
and (O, respectively. 

In summary, to derive the three-body bound state with 
BOA, we should first find the ground-state solution t/j of 
the instantaneous eigen-equation (j6j) of the light atom 
5, and then solve the effective eigen-equation ([7]) of the 
heavy atoms Ai^2 where the instantaneous eigen-energy 
Vef[{zi^Z2) of tjj plays a role as interaction potential be- 
tween Ai and A2. Therefore, the BOA provides a simple 
and clear physical picture for the three-body problem, 
i.e., the light atom B induces an effective interaction be- 
tween the two heavy atoms, which determines the prop- 
erties of the three-body bound state. With this picture, 
one can perform not only quantitative calculations but 
also qualitative discussions for the appearance and fea- 
tures of the trimer states when the potential function 
Kff(^i, ^2) is known from (j6]). This is a major advantage 
of the BOA approach. 

In the end of this subsection we emphasize that, since 
in BOA the transitions between different solutions of the 
instantaneous eigen-equation (|6]) is neglected, this ap- 
proximation can only be used when the gap between 
Kff (^1, ^2) and other eigen-energies of (|6]) [with boundary 
conditions (|4]) and (j5j)] is large enough. In the cases where 
Vef[{zi^Z2) is close to the lower bound of the continuous 
spectrum, the application of BOA may be questionable. 



Effective interaction between the two heavy 
atoms 



In the discussion above we outline the procedure for the 
derivation of the three-body bound states with BOA. In 
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this subsection we solve Eqs. (j4][6]) to calculate the in- 
stantaneous eigen-state ip{rB^ ^i-, ^2) of the light atom 
and the light- atom- induced effective potential leff(^i, ^2) 
between the two heavy atoms. 

A straightforward calculation shows that the lowest 
ground state solution (up to a normalization factor) 
of Eq. ([6j) and the corresponding energy leff(^i,^2) are 
given by 



Kff (2^1,^2) 



p — K{ri2)riB p — i^{ri2)r2B 

+ e(ri2)— ^ ; (8) 



riB 



r2B 



(9) 



where ri2 = y^l + {zi — Z2)^ is the distance between A\ 
and A2. 

Substituting the expression of '0(rB,^i,^2) into the 
Bethe-Peierls boundary conditions (j4]) and (j5j), one can 
derive the values of k and ^ in terms of the distance 
{z\ — 2^2), and then obtain expressions for V^(rB, ^1, 2^2) 
and Kff(2^i, 2^2). Notice that, as a bound state, the wave 
function '0(rB,^i,^2) must approach zero in the limit 
r\B — > 00 or T2B 00. Therefore, the condition > 
must be satisfied when we solve the equations of and 

According to Eq. (|9]), the effective potential Vef[{zi^ Z2) 
is a function of distance 2^12 = 2:1 — 2:2 between the two 
heavy atoms along the axial direction of ID tubes. Then 
the wave function ^(2:1,^2) in the total wave function 
(|3]) is also a function of 2^12, indicating the translational 
symmetry along the 2;- axis. From now on, we rewrite 
Kff (^1,2:2) as 14ff(2:i2), and (/)(^i,^2) as (I){zi2), and write 
Eq. ^ as 



2m* dzi2 



^{Z12)=E^{Z12). (10) 



where = 17111712/ {mi -\-17i2) is the reduced mass of the 
two heavy atoms. From Eq. (p!Q|) . we can see clearly that 
Vef^ serves as an effective interaction between the two 
heavy atoms Ai^2, and determines the existence and be- 
havior of the three-body bound states. Next, we discuss 
the feature of V^ff in different parameter regions. 



ai 
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a > 



In this case the two heavy atoms Ai^2 have the same 
positive scattering length with the light atom B. Since tjj 
is the ground-state solution of Eq. (|6]), a straightforward 
calculation shows that in this symmetric case we have 
^ = 1 and given by the equation 



ri2 



(11) 



This equation can be solved analytically, leading to. 
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FIG. 2: (color online) The regularized effective potential 
U{a, Z12) between the two heavy atoms Ai,2 with scattering 
lengths ai = a2 = a = 0.2 (red solid line with open circles), 1 
(blue solid line) and 00 (black dashed line). We also plot the 
effective potential Vef¥{zi2) for ai = a2 = — 5 (green dashed- 
dotted line). The natural unit h = niB — L = 1 is used 
throughout this paper. 



where W{z) is Lambert W function or the principle root 
of equation z = We^ . Substituting the result ([12]) into 
Eq. (j9j), we finally obtain an analytic expression of the 
effective interaction between the two heavy atoms: 



1 

2^' 



Vei^{zi2) = U (a; 2^12) 
where the regularized part U (a; Z12) is given by 



(13) 



U (a;zi2) 



1 



(14) 
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which approaches zero in the limit |2;i2| 00. Therefore, 
the characters of bound states are essentially determined 
by the behavior of U (a; 2^12). 

With the knowledge of the W function, we can easily 
find that when a > 0, {a; Z12) is a pure symmetric 
potential well with 



U (a, Z12) = U (a, -Z12) < 0. 



(15) 



In Fig. [21 we plot U (a; 2^12) for a set of typical values 
of scattering lengths. It is clearly shown that U (a; 2^12) 
provides a simple ID potential well for the two heavy 
atoms. This behavior guarantees that there exists at least 
one bound-state solution (j) of Eq. ([TQ|) . and then the total 
system has at least one three-body bound state. 

Intuitively speaking, one would expect that the atom- 
atom interaction effect be most significant when the scat- 
tering length takes infinite value. However, we find from 
Fig. 2 that the depth of the effective interaction U (a; Z12) 
takes a maximum value when a = 1 in our natural unit. 
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FIG. 3: (color online) The depth D{a) of the regularized part 
U(a,zi2) of the effective interaction between the two heavy 
atoms in the case ofai=a2 = a>0. Notice that D(a) 
takes a maximum value at a = 1, indicating a new resonance 
behavior for mixed dimensional systems. 



rather than a ^ +oo. This observation suggests that the 
light-atom-induced interaction between the two heavy 
atoms Ai^2 is most significant when the scattering length 
between a single heavy atom and the light one equals to 
the distance separating the two ID tubes. This novel 
property can be considered as a kind of resonance effect 
given by the special configuration of mixed dimensional 
systems. 

This resonance effect can also be proved analytically 
with the character of the W function. For any given 
value of a, the potential U (a, Z12) fias only one minimum 
point, which is localized at the origin Z12 = 0. Thus, the 
depth of the potential well takes the form 



D{a) = -/7(a,0) 



1 



W 



-We 



-i/'^) .(16) 



It is easy to show that D (a) takes the maximum value 
when a = 1. In Fig. [3] we plot the potential depth as 
a function of 1/a, exhibiting the resonance signature at 
a = 1. 



ai 



a2 = a < 



In this case, by substituting Eq. (j8j) into the Bethe- 
Peierls boundary conditions dU and (|5]) , we also get ^ = 1 
and Hi given by Eqs. (pT]) and ([12]) for ri2 < How- 
ever, for ri2 > |a|, there is no positive solution of Eq. 
(pTj) for tv. This suggests that the Schrodinger equation 
([6]) with Bethe-Peierls boundary conditions dU and ([5]) 
do not support any instantaneous bound state t/j of the 
light atom 5, and then one cannot derive any effective 
interaction for the two heavy atoms ^41^2 within BOA. 
As a consequence, when the scattering length \a\ < 1, 
there would be no three-body bound state since the con- 
dition ri2 > |a| is satisfied with arbitrary ID distance Z12 



between the two heavy atoms. 

On the other hand, when \a\ > 1, the BOA can give 
the effective interaction potential 



Kff (2:12) 



W 



12 



(17) 



provided that \zi2\ < \/ d?' — \ or ri2 < In the outer 
region of \z\2\ > V a'^ — 1, the potential takes zero value 
as Vef^ (^12) = 0. In Fig. O we also show ]4ff (^12) with 
negative scattering length. 

We would like to emphasize that BOA can only be used 
when Vef[ is well-separated from the continuous spectrum 
of the Schrodinger equation (|6]). This criteria is actually 
broken in the region ri2 ~ |a| or Z12 ~ Va'^ — 1, where 
we have ]4ff (^12) 0. Then the effective potential is not 
apphcable in these regions. Fortunately, if the potential 
is deep enough, the ground-state wave function (j) of the 
heavy atoms ^41^2 would be mainly localized in the region 
Z12 ~ or ri2 <C |a|, where BOA is applicable. Thus, the 
ground-state wave function and its corresponding binding 
energy obtained from BOA is still reliable. Notice that in 
this negative scattering length regime, ^41^2 and B cannot 
form any two-body bound state, hence the appearance of 
a three-body bound state is a non-trivial universal phe- 
nomenon. 



3. < ai < a2 or a2 < < ai 

Now we consider the general cases where the scattering 
lengths ai and a2 are different. In these cases one can also 
derive the values of ^ and hi by substituting the expression 
(|8|) into the Bethe-Peierls boundary conditions dH and 
(j5j). When < ai < a2 or a2 < < ai, we know that in 
the limit ri2 00, that is the two heavy atoms are far 
away from each other, the instantaneous ground state of 
the light atom B is the two-body bound state of B and 
Ai. Considering the expression (j8j) of the instantaneous 
bound state, we have 

e(ri2^oo)=0. (18) 

With the help of this condition, we obtain the result 



-A+ a/A2 + 4e-2'^'-i2/r2 



where 



A.l-1>0. 

ai a2 

Then the value of is given by 



(19) 



(20) 



12 



1 

ai 



-■ (21) 



By solving Eqs. (jT9|) and (j2T]) numerically, we can obtain 
the values of ^ and and then the effective potential 



6 



14ff- It is easy to show that V^ff < for all values of 2^12. 
Therefore, there is also at least one three-body bound 
state. When a2 < < ai , although the atoms Ai and B 
can form a two-body bound state, there is no two-body 
bound state for A2 and B. In this sense the existence of a 
three-body bound state is also a non-trivial phenomenon. 



4- ai < a2 < 

In this case a straightforward calculation shows that 
the values of ^ and are also determined by Eqs. ([T9j) and 
(|2T]) . Nevertheless, similar to the case of ai = a2 = a < 
0, there are also some regions where the instantaneous 
bound state i/j does not exist. Specifically, we can define 
a critical distance as 



m =3.33 



m =9.5 



12 



2 

ai 



1-1/2 



(22) 



with A defined in ([2Q|) . It is apparent that when ri2 > 
r22, we cannot find any real tv which satisfies Eq. (f2T]) . 
In this sense, can be understood as the range of the 
effective interaction between Ai and A2. When this range 
is smaller than the distance between the two ID tubes, 
i.e. ri2 < 1, the two heavy atoms are always separated 
far enough such that the BOA does not give any effective 
mutual interaction. On the other hand, when rj2 ^ 1 the 
effective potential of the two heavy atoms can be defined 
as 



eff 



72; 
0; 



1 < ^12 < r 
ri2 > rl2- 



12 



(23) 
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This potential is also not reliable in the region ri2 ^ r 
where the condition for BOA is broken. However, as 
shown blow, this approach could lead to the existence 
of a bound state wave function (j) which takes negligible 
value in this questionable region, such that the discussion 
within BOA remains valid. Since the negative scattering 
lengths do not support any two-body bound states, the 
existence of such a three-body bound state in this region 
is of great interest. 



III. THREE-BODY UNIVERSAL BOUND 
STATES IN 1D-1D-3D SYSTEMS 

In the previous section, we have obtained the instanta- 
neous bound-state wave function i/j of the light atom B 
and the effective interaction potential leff(^i2) between 
the two heavy atoms. We have shown that Vef[ is most 
significant when the two-body scattering length is reso- 
nant with the distance between the two ID tubes. In this 
section we derive the wave functions and binding ener- 
gies of the relevant three-body bound states, and further 
confirm the observation of this new resonance effect. 





FIG. 4: (color online) The binding energy of the ground 
three-body bound state in the 1D-1D-3D system with reduced 
mass of the heavy atoms m* = 3.33 and 9.5. These values 
correspond to the cases of (Ai = A2 =^°K, B =^Li) and 
{Ai = A2 =^^^Cs, B ='^Li), respectively. 




FIG. 5: (color online) The binding energy Esh of the ground 
trimer state as a function of 1/a with ai = 02 = a. The 
reduced masses used in this plot are m* = 3.33 (black solid 
line), 9.5 (blue dashed line) and 00 (red solid line with open 
circles) , respectively. 



In Fig. m the binding energy Es\) of the ground trimer 
state is plotted as a function of 1/ai and l/a2 with heavy- 
atom reduced masses = 3.33 and 9.5 in the natural 
unit. These values correspond to the cases of {Ai = 
A2 ='^°K, B =^Li) and {Ai = A2 =^^^Cs, B =^Li), 
respectively. Here, the binding energy £^3b is defined as 
the energy gap between the three-body ground state E 
and the threshold of the effective interaction, i.e. 



^3b = Kff(00)-^. 



(24) 



From Fig. HJ we notice that a three-body bound state 
exists for a wide range of positive and negative scattering 
length combinations, as discussed in the previous section. 
Nevertheless, the binding energy reaches a peak value 
when the two scattering lengths ai and a2 are close with 
each other, especially in the region around ai ~ a2 ~ 1. 
This observation is consistent with the expectation out- 
lined in the previous section, which shows that when 
ai = a2 = a, the effective potential well for A1-A2 inter- 
action is deepest as the scattering lengths are resonant 
with the distance between the two ID tubes a = 1. 
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FIG. 6: (color online) The binding energy of the ground 
trimer state as a function of the reduced mass m* with ai = 
a2 = a — 1 (blue solid line with open triangle) and ai = a2 = 
a = oo. 



To further investigate the relationship between the 
binding energy and the two-body scattering lengths, we 
focus on the case of ai = a2 = a, and illustrate in Fig. [5] 
the binding energy in terms of 1/a with respect to differ- 
ent reduced masses of the two heavy atoms. One sig- 
nificant feature of this result is that the resonant behavior 
is present for all different reduced masses, i.e. the binding 
energy of the ground trimer state reaches its maximum in 
the region around a = 1. Besides, we also notice that for 
a given two-body scattering length, the binding energy 
increases with reduced mass m*, and approaches to an 
asymptotic value in the limit oo. This tendency is 

also confirmed by Fig. [6] where the binding energies for 
ai = a2 = a = 1 and ai = a2 = a = oo are plotted as 
functions of the reduced mass m*. 

The three-body bound states in the 1D-1D-3D systems 
with ai = a2 are also discussed in Ref. j36| within an 
effective field theory or the exact solution of three-body 
Schrodinger equation. In Fig. [71 we compare our BOA 
results of the ground trimer state energy with the exact 
expression given by Ref. [SG^ for = 3.33 and ai = 
a2 = a. Notice that the BOA results are very close to the 
exact solution around the resonance point a = 1 for such 
a rather small mass ratio. This consistency suggests that 
the BOA approach is reliable provided that the three- 
body bound state energy is away from the threshold. 



FIG. 7: (color online) The binding energy of the ground three- 
body bound state in the 1D-1D-3D system with reduced mass 
m* = 3.33 and scattering lengths ai = a2 = a. Here, we plot 
the results given by the BOA (red solid line with open circles) 
and by the exact solution of the Schrodinger equation [3y] 
(blue solid line). Notice that the BOA can give reliable results 
provided that the binding energy of the trimer state is away 
from the threshold. 



When the masses of ^41^2 is much larger than that of 
5, the system can also be treated via BOA. The wave 
function ^ of the possible three-body bound state also 
takes the factorized form as in Eq. ([3]), i.e., ^ = 0?/; 
with tjj the instantaneous bound state of the light atom 
B. In this case, the instantaneous energy of tjj serves as 
an effective 2D interaction between the two heavy atoms, 
and can be obtained by re placing the argument Zi — Z2 
in Feff(^i - Z2) with p = ^/{yi - ^2)^ + {zi - Z2Y . Fol- 
lowing the same procedure as outlined in Sec. II, we can 
show that in the case of ai = a2 = a, the depth of the 
2D effective potential also takes its maximal value when 
a = 1 in the natural unit. This observation indicates that 
the resonance phenomenon also exists in the 2D-2D-3D 
configuration. 

Notice that the 2D-2D-3D geometry is invariable un- 
der a rotation along the x-axis. This SO (2) symmetry 
thus leads to the conservation of the x-component angu- 
lar momentum of A1-A2 relative motion. Therefore, the 
wave function (j) in the three-body bound state ^ can be 
expressed as 



E 



(25) 



IV. THREE-BODY UNIVERSAL BOUND 
STATES IN 2D-2D-3D SYSTEMS 

The discussion on 1D-1D-3D systems outlined in the 
previous section can be directly generalized to other 
mixed-dimensional configurations. In this section we con- 
sider a 2D-2D-3D system [Fig. 1(b)] where the two heavy 
atoms Ai^2 are trapped individually in two 2D confine- 
ments, localized in the planes of x = ±L/2. The light 
atom B is also assumed to move freely in the 3D space, 
we also adopt the natural units with h = tub = L = 1. 



where tan6> = {zi — Z2)/{yi — ^2) is the polar angle of 
Ai^2 relative motion in the y-z plane, and the radial wave 
function (/)^(p) satisfies the 2D Schrodinger equation 



1 



2m* 



d^^-pd-p-^)^""^^^'^ 



(26) 



Here, the quantum number £ = 0, ±1, ±2, ... indicates the 
relative angular momentum of Ai^2 along the x-direction. 
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FIG. 8: (color online) The binding energy of the ground three- 
body bound state in the 2D-2D-3D geometry with reduced 
mass of the heavy atoms m* = 3.33 and 9.5. 




1/a 



FIG. 9: (color online) The binding energy of the ground three- 
body bound state in 2D-2D-3D geometry with reduced mass 
m* = 3.33 and scattering lengths ai = a2 = a. Here, we 
plot the results given by the BOA (red solid line with open 
circles) and by an effective field theory [s^] (blue solid line), 
and find good agreement provided that the binding energy is 
away from the threshold. 



The pure ground state of the system occurs in the channel 
£ = 0. 

The radial equation (|26|) can be solved numerically as 
in the 1D-1D-3D case. For the ground zero-angular mo- 
mentum channel ^ = 0, we also find three-body bound 
states with reduced mass = 3.33 and 9.5. The binding 
energy of the ground trimer state is illustrated in Fig. [8]in 
terms of 1/ai and l/a2. Notice that the binding energy is 
significantly amplified in the parameter region ai ~ a2, 
and reaches its maximum when the scattering lengths 
are resonant with the 2D surfaces spacing ai ^ a2 ^ I. 
Besides, the binding energy also increases with reduced 
mass of the two heavy atoms. In Fig. [9] we also com- 
pare the BOA results with the exact expression [32| for 
the case of = 3.33 and ai = a2 = a, and find good 
agreement when the trimer binding energy is away from 
the threshold. All these features are analogous with the 
case of 1D-1D-3D geometry. 



1D-1D-1D-3D system 




X 

-L/2 L/2 



FIG. 10: (color online) The 1D-1D-1D-3D system with three 
heavy atoms Ax^A^ and A^ confined in three ID tubes and 
the light atom B moving freely in 3D. 



V. FOUR-BODY UNIVERSAL BOUND STATES 
IN 1D-1D-1D-3D SYSTEMS 

From the discussion in the previous sections, we no- 
tice that the BOA works well throughout a wide range 
of scattering length for a fairly small mass ratio of about 
6, provided that the binding energy of the bound trimer 
state is away from the threshold. This observation sug- 
gests that this approach can be directly applied to mixed 
dimensional systems with more than three atoms, and 
to give reliable results for few-body bound state energy 
when it is sizable. In this section, we consider as an ex- 
ample the 1D-1D-1D-3D system with three heavy atoms 
Ai, ^2, and A^ trapped individually in parallel ID tubes 
and a single light atom B moving freely in 3D. 

We consider the configuration of three ID tubes ar- 
ranged along the z direction, and intersect with the x-y 
plane at [x = ±L/2,y = 0) and {x = xo^y = ^o), as 
shown schematically in Fig. |Vl The three intersection 
points form a triangle in the x-y plane. Since the system 
properties are invariant under different length scales, we 
assume that L is the shortest side of the triangle, and use 
it as the length unit L = 1 in the following discussion. 

The quantum states of such a system can be described 
by the wave function ^(r^; ^i, ^2, ^3), where Zi is the z- 
coordinate of the heavy atom Ai , and tb is the coordinate 
of the light atom B. Within the BOA, the wave function 
^ can be separated as 

^(fe; zi,Z2, Z3) = ^(^1, ^2, ^3)^(^8; zi,Z2, Z3). (27) 

Here, is the wave function of the instantaneous bound 
state of the light atom, which is given by the Schrodinger 



equation 



with the Bethe-Peierls boundary conditions 

^ 0) oc f 1 - ^) + OinB). 



(28) 



(29) 



Here, tib and are the distance and mixed-dimensional 
scattering length between the atoms Ai and 5, respec- 
tively. 

The ground state of ^/^ can be obtained by solving the 
eigen-equation ([28|) for a give set of (^1,2:2,^3), which 
takes the form 




FIG. 11: (color online) The regularized effective potential 
[/(a; X, y) for two body scattering length ai = a2 = as = 
a = 1 in the 1D-1D-1D-3D system with equilateral triangle 
configuration. 



^(fB\Zx,Z2,Zz) 



riB 
^2 



C2- 



r2B 



C2- 



r2B 



Vef^{zi,Z2,Z3) —. 



(30) 



The parameter is determined by the boundary condi- 
tions ([29]) , leading to 



ri2 



C3- 



C2t^ - C3- 



ri2 

C2- 



ri3 



r23 



1 

1 

ai 
C2 . 

a2' 



ri3 



r23 



C3 

C3K. = — . 

«3 



(31) 



A numerical solution of these equations hence gives the 
effective potential Vef^ among the three heavy atoms. By 
imposing the potential to the Schrodinger equation 



E 

i=l,2,3 



1 

2mi dzf 



Vei^{zi,Z2,Z3) 



0(^1,^2,^3) 



= E^{zi,Z2,Z3), (32) 



we can obtain the energy E for the four-body bound 
states. 

From now on, we focus on the special case of ai = a2 = 
a3 = a and mi = m2 = ms = m, that is the scattering 
lengths and the masses of the three heavy atoms are all 
the same. This is also the most relevant case for experi- 
ments, where atoms trapped in low dimensional traps are 
of the same species. Since the system is translationally 
invariant along the 2;-direction, we define a new set of 
coordinates 



X 

Y = Z3 



Zl - Z2, 

Zl + Z2 



(33) 



and calculate the effective potential Vei^{X^Y) in these 
new variables 



Veif{X,Y) = U{a;X,Y) 



1 

2^' 



(34) 



where U{a]X^Y) is the regularized part. 

We first consider the special geometry where the three 
ID tubes are arranged equidistantly to form an equilat- 
eral triangle in the x-y plane (i.e., xq = and yo = 
V3/2). In Fig. [m we show the regularized effective po- 
tential /7(a;X, y) for scattering length a = 1. Notice 
that the effective potential acquires its global minimum 
at (X = 0,7" = 0) or zi = Z2 = ^3, that is the three 
atoms are staying in a surface perpendicular to the ID 
tubes and forming an equilateral triangle. Besides, we 
also observe three energy potential valleys, which corre- 
spond to the cases where the distance between two of the 
three atoms equals to 1. 

The same phenomenon can also be observed for other 
values of scattering length a ^ 1. In fact, the effective 
potential U{a;X^ Y) always reaches its minimum at {X = 
0,F = 0). However, the potential is deepest only when 
the scattering length a = 1. In Fig. [121 we show the 
depth of the effective potential well as a function of a, 
which reaches its maximum at a = 1. This result suggests 
that the resonance we observed in three-body problems 
as discussed above also occurs in the four-body system. 

With the knowledge of the effective potential, we can 
numerically solve the Schrodinger equation ([32]) to obtain 
the eigenenergies of four-body bound states. In the new 
set of variables X and F, this equation can be rewritten 
as 



1 



3 



I dX^ 4m dY^ 



Kff(x,r) 



(/)(x,r) 



E^{X,Y), 



(35) 



where (/)(X, Y) is the wave function of the heavy atoms. 
As in the three-body calculation, the binding energy of 
the tetramer states is defined as the difference between 
the eigenenergy E and the effective potential energy for 
X ^ 00 and y ^ 00, 



E4h = V;ff(oo,oo) - E. 



(36) 



The binding energy of the ground four-body bound 
state for different values of a is plotted in Fig. [131 where 
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FIG. 12: (color online) The depth D{a) of the regularized 
part U(a; X,Y) of the effective interaction. In this plot, we 
consider the case of ai = a2 = as = a in the 1D-1D-1D-3D 
system with equilateral triangle configuration. Notice that 
D{a) takes maximum value at the resonance point of a = 1. 
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FIG. 13: (color online) The binding energy E^h of the ground 
four-body bound state in the 1D-1D-1D-3D system with equi- 
lateral triangle configuration, for different values of a with 
reduced mass m* = 9.5 (blue solid line with circles) and 3.33 
(green solid line with triangle). 

we consider two mass ratios as in the previous discus- 
sion. Notice that the binding energy reaches its max- 
imum near a = 1, as we expected from the effective 
potential. This result confirms the appearance of the 
resonance phenomenon in the four-body system. 

Up to now, we consider only a special configuration of 
1D-1D-1D-3D geometry where the three ID tubes form 
an equilateral triangle, and observe a resonance phenom- 
ena for tetramer binding energy as the scattering length 
gets close to the mutual distance between ID tubes. An 
intuitive expectation is that this most symmetric config- 
uration should be the case of maximal resonance, for the 
scattering length can be resonant with any two of the 
three atoms. In order to demonstrate this idea, we con- 
sider general configurations of the three ID tubes, such 
that they form a triangle of arbitrary shape with three 
sides L = 1, Li and L2 (see Fig. 10). For the sys- 
tem properties are invariant as scaled with length, we 
assume L = 1 to be the shortest side of the triangle. We 
further take the scattering lengths ai = a2 = as = 1. 
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FIG. 14: (color online) The depth of the effective potential 
U{a; X, Y) in the 1D-1D-1D-3D system with ai = a2 = as = 
1 and the inter-tube distances Li — 1 and Li,2 defined in Fig. 
10. 

In Fig. [m we show the depth of the effective poten- 
tial [/(a; X, Y) for arbitrary arrangement of the three ID 
tubes. It is clearly shown that, the depth of the effective 
potential takes its maximum value when Li = L2 = 1 or 
the ID tubes have the configuration of equilateral trian- 
gle. This is consistent with our expectation that maximal 
resonance appears in this most symmetric configuration. 



VI. BOA FOR MANY-BODY PROBLEMS IN 
MIXED-DIMENSIONAL SYSTEMS 

In the previous sections, we study the three-body 
and four-body bound states in mixed-dimensional sys- 
tems within the BOA. Now we generalize this approach 
to mixed-dimensional problems with arbitrary heavy 
atoms trapped individually in ID or 2D confinements, 
while a single light atoms moving freely in the 3D space. 
In such a configuration, the wave function of the possible 
few-body bound states takes the form 

^I^{rB;s)=^{s)^{rB;s), (37) 

where s = (ri,r2, ..,r/v) are the ID or 2D coordinates of 
the heavy atoms Ai, ^2, Aat, and fs is the coordinate 
of the light atom B. As in the previous sections, tjj{fB] s) 
is the wave function of the instantaneous bound state of 
the light atom, which is determined by the Schrodinger 
equation 

- \^IWb-. s) = V,^{s)^{rB-, s) (38) 
with Bethe-Peierls boundary conditions 

^{ns ^ 0) (X f 1 - ^\ + 0(r,B). (39) 

Here, Vis is the distance between the atoms Ai and B. 
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By solving Eq. (|38|) . we obtain the general form of the 
instantaneous bound state 



riB 



N 

E 

i=2 



Ci- 



(40) 



where the value of and the coefficients Ci are given by 
the equations 



1 ^ 

ci e 
— = i^ci 



Ci- 



TiB 



riB 



N 

E ' 



riB 



(41) 



(42) 



From the equations above, we can solve for the value of 
Hi in terms of the coordinate s of the heavy atoms, and 
then obtain the instantaneous wave function of '0(rB; s) 
and the effective interaction among the heavy atoms 



"T 



(43) 



Finally, the heavy- atoms wave function (j){s) of the few- 
body bound state is given by 



^ 1 



s) = E(t){s) (44) 



with rui the mass of the heavy atom Ai. 



VII. CONCLUSION 

In this manuscript we show our BOA-based results 
on the stable three-body or four-body bound states in 
mixed dimensional systems with N > 2 heavy atoms in- 
dividually trapped in different ID or 2D confinements, 
while a single light atom moving freely in the 3D space. 
The BOA approach can provide a clear physical pic- 
ture with a well-defined effective interaction among the 
heavy atoms. We show that in mixed dimensions, the 
three-body or four-body bound states can occur within 
a broad range of two-body scattering lengths, as the Efi- 
mov states in 3D. Nevertheless, the binding energy of the 
ground bound state reaches its maximum value when the 
two-body scattering length gets close to the distance be- 
tween the low-dimensional traps. This is due to a new 
resonance phenomenon in mixed dimensions, where the 
effective interaction among the heavy atoms acquires a 
deepest potential well under the resonant condition. The 
feasibility of this BOA approach is then confirmed by a 
direct comparison with exact results in 1D-1D-3D and 
2D-2D-3D configurations, hence suggests a possible ex- 
tension in the problems with more than three atoms in 
mixed dimensions. 
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Appendix A: The Bethe-Peierls Boundary Condition 
for BOA in Mixed-Dimensional Systems 

In this appendix, we derive the Bethe-Peierls boundary 
condition [e.g., Eqs. (g]), (|5|), and (^] used in the 
Born-Oppenheimer approach for the mixed dimensional 
systems. For simplicity, here we consider the system with 
one heavy atom A confined in a ID trap which is arranged 
along the z-axis, plus a light atom B moving freely in 3D. 
The generalization to other cases is straightforward. 

The expression of Bethe-Peierls boundary condition 
should be derived from the asymptotic behavior of the 
two-body wave functions. As a first-principle discus- 
sion, we first take into account the 3D motions of both 
atoms A and 5, and then reduce our result in the mixed- 
dimensional model where only the motion along the z di- 
rection is considered for atom A. The total Hamiltonian 
of the two atoms is given by 

HAB=TAz^TA±^VA±^TB^VAB{rAB)^ (Al) 

Here, the kinetic energy of atom A along the z direction 
is given by 



Taz = 



1 

2mA dz\ 



(A2) 



with tua the mass of atom A, and fi=A,B = (^i, Vi^ ^i) the 
coordinate of the corresponding atoms. The transverse 
kinetic energy Ta± of atom A and the total kinetic energy 
Tb of atom B are defined as 



Ta± = 
Tb = 



1 



2mA \dx\ dy\ 

-hi. 



(A3) 
(A4) 



Here we use the natural unit h = mB = I- In the Hamil- 
tonian ()Aip we also have the transverse harmonic poten- 
tial 



Va± 



(A5) 



with frequency uj±, and the atom-atom interaction po- 
tential Vab {^ab) which is a function of the distance be- 
tween the two particles tab = \^a — ^bI- We further 
denote the effective range of the interaction potential as 
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, such that we have Vab {i^ab ) ~ in the region of 
Tab > r^. 

When the confinement of the heavy atom A is strong, 
the transverse motion of atom A in the x-y plane is much 
more rapid than its motion along the z direction. There- 
fore, we need to consider both the position re of the 
light atom B and the transverse coordinates (xa^Va) of 
the heavy atom A as fast degrees of freedom. Only the 
longitudinal coordinate za of atom A is treated as the 
slow variable. 

Within the BOA, the total wave function of the system 
takes the form 

^ (rA, vb) = (j) {za) ij{fB,XA, Va; za), (A6) 
where ^{^{rB^XA, z/a'^za) is given by the eigen-equation 

Hp {za) i^{rB,XA, Va; za) = E {za) ilj{fB,XA, Va; za) 

(A7) 

of the Hamiltonian 

Hf {za) = Ta± + Va± ^Tb^ Vab {tab) (A8) 



with fixed values of za- To solve Eq. (|A7|) . we expand 
the solution tjj with eigen-states of the transverse Hamil- 
tonian Ta± + Va± of atom A 

oo 

tlj{fB,XA,yA;ZA) = X^0n V^n {tB] Za) • (A9) 

n=0 

Here, (j)n {xa^Va) is the n*^ eigen-state of Ta± + Va±' 
Considering the translational symmetry along the z-axis, 
we take za = 0, and the relevant wave function ipn {tb] 0) 
of the light atom B is given by 



[TB + (n + l)cj^]V^n + 



^Vnm {vb)"^ 
m 

= E{0)^n^ (AlO) 



Here, the matrix element of the interaction potential 
takes the form 



Vnm {tb) 



j dxAdyA(t>n {^a.Va) 

>^VAB{rAB)(t>m{xA,yA) ' (All) 



Therefore, the eigen-equation (jA7p or (|AlQp can be 
solved via a multi-channel scattering theory of atom 5, 
with the transverse states (pn {tb] za) of atom A serving 
as the scattering channels. In the low-energy case with 
uo^ < E < 2a;^, the ground channel with the transverse 
state {xa^Ha) assumes the only open channel. 

Now we consider the asymptotic behavior of the wave 
function in the long-distance limit with \fB\ ^ 
where l±_ = y^l/ {mA^±) denotes the characteristic 
length of the transverse confinement. In this region, the 
mutual distance tab between the two atoms would be 
much larger than the effective range of the interaction. 



such that we can neglect the term Vab in Eq. ()A8p . Ac- 
cording to the scattering theory, in such a region the wave 
function tjjn (^b;0) in the close channels with n > de- 
cays exponentially with |rB|, and can be safely neglected. 
The wave function i/jq {vb^O) in the open channel takes 
the form 



1=0 m=-l 



[ji {k Irsl) + kfi,^ [k) h\+^ {k Ifsl)) , (A12) 



where k = ^/2{E — oj±)^ Yim (j)) are the spherical har- 
monic functions of the azimuth angles {9b-,4>b) of r^, 
ji {z) is the Riccati-Bessel function, and h['^^ {z) is the 
Riccati-Hankel function. The coefficients Ci^rn are given 
by the boundary condition, while the scattering ampli- 
tudes fi^rn {k) are determined by the effective potential 
Vnm {^b) defined in ()Alip . In the low-energy case with 
small /c, we can neglect all the high-partial wave scat- 
tering amplitudes fi^m{k) with / > 0, and approximate 
the 5- wave scattering amplitude /o,o {k) with /o,o {k = 0). 
Then the long-distance behavior of wave function ip be- 
comes 



ilj{fB,XA,yA;0) 00 {xA,yA)'ipo {rB;0) 



00 {xa.Va) 



f sin {k 



I 



k \ fB 



yBj(PB 



klTBl 



re I) - kaABe 



Ji {k\rB\) 



ik\rB\ 



I— I rn= — l 

with the scattering length qab defined as 
aAB -/o,o {k = 0). 



(A13) 



(A14) 



The expression (|AT3)) impl ies that in the "intermedi- 
ate" region of 



[r*,l±] < Irel < 



(A15) 



the behavior of tjj takes the form of 



V^(rB,^A,yA;0) - 00 {xa^Va) (^l - . (A16) 

Therefore, we can replace the real interaction potential 
Vab {'^ab) in (|Aip with a Bethe-Peierls-type boundary 
condition 



lim V^(fB,XA,^A;0) oc 00 ( 1 - 7^ 



(A17) 



Under this boundary condition, the solution of the eigen- 
equation 

[Ta± + Va± + Tb] Wb, xa, y^; 0) 

= E'tl^{rB,XA,yA;^) (A18) 
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takes the form of Eq. ()A13p for all |rB | 7^ 0, and becomes 
a reasonable approximation for the solution of (|A7p . 

In this reduced mixed-dimensional model, the trans- 
verse coordinates (xa^Va) of the heavy atom A is taken 
to be fixed values of (0, 0). Together with the assumption 
= 0, we have 

\fB\=rAB. (A19) 

and then the boundary condition (jA17p can be expressed 
as 



with 



lim (re; 0) oc (1 

VAB^O 



TAB 



(A20) 



Here, i/j {vb^O) is the wave function of the light atom B 
with the position of atom A fixed dit za = 0. For non-zero 
za, the condition (|A2Qp can be generalized to 



lim ijj (re; 2;^) oc ( 1 

VAB^O 



O^AB 
TAB 



(A21) 



That is the Bethe-Peierls boundary condition used in the 
BOA discussed in the main text of this manuscript. 

We notice that there is another type of Bethe- 
Peierls boundary condition as discussed in Ref. ^| , 
where the total wave function ^ of the reduced mixed- 
dimensional two-body problem is assumed to satisfy the 
condition 



lim ^ oc 

Dab^O 



1 



Dab J 



(A22) 



Dab 



ruA 



This condition is slightly different from our result of Eq. 
()A2ip . The difference can be understood by noticing 
that when solving for the wave function of atom B under 
BOA, we fix the position of the heavy atom A, such that 
the relevant Bethe-Peierls boundary condition (jA22p be- 
comes isotropic. It is pointed out that, in the limit of 
rriA ^ 1, the condition ()A22p approaches to ()A2ip and 
we have qab = <^eff- Therefore, we approximate qab as 
tteff when comparing the BOA results with the effective 
field theory [36[. 

It is straightforward to generalize the discussion above 
to more general cases with heavy atoms Ai^...^An 
individually confined in N low-dimensional traps, and 
one light atom B moving freely in 3D. In that case, we 
can fix the positions of the heavy atoms under BOA, and 
use the Bethe-Peierls boundary condition 



lim (fAjTB) (X ( 1 - — 



(A24) 



to solve the Schrodinger equation of the light atom. Here, 
riB is the distance between the heavy atom Ai and the 
light atom B. That is the approach we used in our main 
text. 
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